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We review some recent developments in the studies of direct CP violation and final state inter- 
actions in weak decays of heavy quarks. 
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I. INTRODUCTION AND MOTIVATION 

Heavy quark decays serve as a powerful tool for testing the Standard Model and, since it probes the quarks of all 
three generations, provide invaluable possibilities to study CP violation. However, the interpretation of experimental 
observables in terms of fundamental parameters is often less than clear. Rare hadronic decays of B mesons, for 
example, proceed through both tree level Cabibbo-suppressed amplitudes and through one loop penguin amplitudes. 
On the one hand, this situation allows direct CP violating effects. On the other, these contributions complicate the 
extraction of Cabibbo-Kobayashi-Maskawa (CKM) angles and, in particular, the angle 7 = arg [—V u dV* b /V c dV* b ] . 

The decay of a heavy hadron produces quarks in the final state. Because of their strong QCD interactions they 
continue to interact after the weak transition took place. Even after they have formed hadrons, there are still strong 
forces between them, and therefore, the problem of the final state rescatterings or final state interactions (FSI) is 
an important part of the physics of nonleptonic 6-decays. The most pronounced effect of FSI is clearly in direct CP- 
violation where one compares the rates of a £>-meson decay with the charged conjugated process. The corresponding 
\Q \ asymmetries between the two decays depend on both a weak (CKM) and a strong rescattering phase provided 
by the FSI. Furthermore, a non-vanishing asymmetery requires two different final states produced by different weak 
■ amplitudes which can go into each other by a strong interaction rescattering. Thus, FSI directly affect the asymmetries 
and their size can be interpreted in terms of fundamental parameters only if these FSI phases are calculable. 

As an example, let us first look at the recently measured combined branching ratios for i? ± — > 7r ± AT and Bd — > 
tt t K ± jjj. In the Standard Model, these decays are mediated by the AB = 1 Hamiltonian, which takes the form 

Qlh' f-i / 2 6 10 \ / 2 6 10 \ 

T r [ v cb v; s £ aQr + ]T CiQi + ]T dot + v ub v u * s ]T dQf + £ c t Qi + £ c&t 1 + h.c, (i) 
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The flavor structures of the current-current, QCD penguin, and electroweak penguin (EWP) operators are, respec- 
. £h tively, Q\ s 2 ~ sqqb, Qf 6 ~ sb^q'q', and 10 ~ sb^2 e q > q' q' , the sum is over light quark flavors. The Wil- 
son coefficients Cj are renormalization scale dependent, 61,2(1 GeV) — 0(1), (73,. .,8,9(1 GeV) = 0(1O~ 2 ), and 
07,8.10(1 GeV) < 0(1O~ 3 ). Let us expand the decay amplitudes of interest according to their dependence on the 
elements of the CKM matrix, 

A(B + -> 7T+A- ) = A+ AS.e*e*+, A(B~ - n~ K°) = A+ A+ s e'^e tS \ 
A(B° - 7T-K+) = A° cs - A° us e^e* s °, A(B° n + K~) = A Q CS - A^e^e" , (2) 

where Sq and 5+ are CP-conserving phases induced by the strong interaction. The first and second terms in each 
amplitude correspond to matrix elements of the first and second terms in 7i e tf (or their Hermitian conjugates), 
respectively. Note that each term is by itself scheme and renormalization scale independent. 

One can asses the expected relative contributions of the operators in H c s to a given exclusive decay mode. The 
electroweak penguin operators are commonly neglected, since the contributions with a sizable Wilson coefficient, 
CqQq, are color suppressed or require rescattering from intermediate states. In this case isospin symmetry of the strong 
interactions leads to the simplification A^ s = A+ . It is now believed that the current-current operator contributions to 
A®' s + are roughly of same order as the QCD penguin operator contributions. The contribution of the current-current 
operators to A^ s is also expected to be of the same order, despite the CKM suppression, because of the large value of 
C2, namely, V u bV* s C2 ~ V c bV* s C3 ..6- However, since for S ± — > the relevant quark transition is 6 — ^ aids, one 



1 



might expect the size of A+ s relative to to be highly suppressed by the small ratio \V u bV* s /V c bV* s \ ~ 0.02. This 
would hold equally for the current-current and penguin operators. If, indeed, r + = A+ s /A+ s ~ \VubV* s /V c bV* s | is a 
good approximation, then there are two important consequences: 

(i) Direct CP violation could be observed, in principle, through the CP asymmetry A^p = Aq P (B + — ► tt + K°), 

dir _ BR(B+ -> tt+K°) - BR{B~ -> tC_ K°) _ 2r+ sin 7 sin 5+ 

CP ~ BR(B+ ^ n+K°) + BR(B- n'R ) ~ 1 - 2r+ cos 7 cos 6+ + r\ ' ( ' 

However, it would be small, A^%{B + -> tt+K^) < C(A 2 ), where A ~ 0.22 is the Wolfenstein parameter. "Hard" FSI 
estimates, where the u quarks in Q\\ are treated as a perturbative loop, give A^ P ~ 1%. 

(n) Model-independent bounds could be obtained for the angle 7 using only the combined branching ratios 
BR(B ± -> 7 r ± K) and 5i?(-B d tt^K^ One can construct the ratio 

„ BR(B° ^ 7T-K+) + BR(B° ^ 7T+K-) ( A° cs \ 2 1 - 2r cos 7 cos 6 a + r 2 

R = TTTTTTTI .„n, . Tr^T = ( TT ) 1 o„ x — ~ > ( 4 ) 



BR{B+ ^tt+K q )+BR(B- ^tt-K°) \A+J 1 - 2r + cos 7 cos <5+ + r 



where tq = A^ s /A° cs . If j4+ s and EWP operator contributions are negligible, the ratio (^) takes the simple form 

R = 1 — 2r cos 7 cos 5 + r o ■ (5) 

The observable R may be minimized with respect to the parameter ro, which (as cos 2 So < 1) leads to the bound 

sin 2 7 < R. (6) 

If true, a stringent bound on 7 would be obtained if the reported by CLEO value i? e xp = 1-00 ± 0.40 Q turns out to 
be smaller than unity within experimental errors 

Rare B decays, like B — > 7rX, are suppressed in the Standard Model by either CKM matrix elements or small 
Wilson coefficients. Thus, these decays are potentially sensitive to New Physics. In the presence of New Physics, a 
large CP asymmetry can be induced, Aq P » 1% and R can be modified to violate the bound (||). The analysis leading 
to (|^), however, explicitly assumes that the CKM angle 7 does not enter the theoretical expression for the charged 
decay amplitudes (|^), i.e. the absence of large contributions from the operators Qi%. This assumption is based on 
the observation that the quark level decay b — > dds is not mediated direcly by Q\\. However, this treatment of the 
dynamics ignores the effects of soft rescattering effects at long distances, which can include the exchange of global 
quantum numbers such as charge and strangeness. In the absence of an argument that parton-hadron duality should 
hold in exclusive processes involving pions and kaons, one must conclude that the long distance physics of meson 
rescattering is not probed by the analysis of the final state rescattering based on perturbative QCD. In addition, 
the rescattering in question is inelastic, despite its quasi-elastic kinematics, and cannot be studied adequately in 
any model of purely elastic final state phases. Let us proceed with our example and investigate the impact of final 
state rescattering on the CP asymmetry Aq P (B ± — > tt^K) and the ratio R. The rescattering process involves 
an intermediate on-shell state X, such that B — > X — > Kix. In particular, we assume that there exists a generic 
(multibody) state Kmr. The charged and neutral channel amplitudes can then be written as 

A{B+ - Krrn) = A n + - Age^e** , (7) 
A(B° -> Kmr) = Af s ~ A£°e i7 e i<5 °" . 

Rescattering contributions, again decomposed according to their dependence on CKM factors, are given by 

A(B+ -» Kmr -» n+K°) = S^A n c + - S^A^fe 1 ^ , 

A(B° -» Kmr -» ?T K+) = S%Af s - S^Af s e^ , (8) 

where S™ is the complex amplitude for rescattering from a given multibody final state to the channel of interest. 
In the limit of isospin symmetry A+ s = A® s , and this equality is not spoiled by rescattering effects. The i = 1,3,4 
rescattering amplitudes can be absorbed into the unknown amplitudes in Eq. (0). 
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Let us assume that the rescatterings of transitions mediated by Qi% are significant enough to dominate A+ s , so 
A+ s e lS+ = J2 n ' an d define e = A+ s /A+ s . Let us also assume that rescattering effects do not dominate the 

overall decay, so we may retain just terms linear in e. Then, A^p of Eq. (||) and R of Eq. take the form 

^ dir _ 2e sin 7 sin 8+ R _ 1 - 2r cos 7 cos S + rg 
CP 1 — 2e cos 7 cos <5+ ' 1 — 2e cos 7 cos 8+ 

Once again, we may extremize R with respect to the unknown ro, 

1- ca* ycc* 8 
1 — 2e cos 7 cos 0+ 

Using the same arguments as before with respect to the strong phases <5o and 5+, we find the new bound 

sin 2 7<i?(l + 2e\/l- J R), or |cos 7 | >y/l-R-eR. (11) 

It is clear that even a small rescattering amplitude e ~ 0.1 could induce a significant shift in the bound on 7 deduced 
from R. For example, the fractional correction to the bound on |cos7| is A = eR/Vl - R. The value of A is a 
strong function of the experimentally observed i? xp, A ~ e for R cxp = 0.65 and A ~ 2e for i? eX p = 0.80 The bound 
deteriorates quickly as R cxp — * 1. In addition, Eq. (^|) implies that similar effect could in principle generate an 
0(10%) CP asymmetry which is significantly larger than the bound Aq P (B ± — > n K) ~ 1%. Therefore, in order to 
understand whether a large CP asymmetry signals New Physics, and whether it is possible to obtain a bound on 7, 
it is imperative to study FSI to obtain an order of magnitude estimate of the effect. 



II. FINAL STATE RESCATTERING 

Final state interactions arise as a consequence of the unitarity of the S- matrix, S^S — 1, and involve the rescattering 
of physical particles in the final state. The T-matrix, defined by S = 1 + iT , obeys the optical theorem: 

ViscT B ^ f = I [(f\T\B) (f\Ti\B)] =l"£{f\T^m\T\B) , (12) 

i 

where Disc denotes essentially the imaginary part. Using CPT in the form (f\T\B)* = {B\T*\f) = (f\T*\B) this 
can be tranformed into the more intuititve form 

(/lT|£)* = £</|<St|i)<i|T|B) . (13) 

i 

Here, the states \i) represent all possible final states (including |/) itself) which can be reached from the state \B) 
by the weak transition matrix T. The right hand side of Eq. ( [T3| ) can then be viewed as a weak decay of \B) into 
|i) followed by a strong rescattering of |i) into |/). Thus, we identify (f\S^\i) as a CP-conserving FSI rescattering of 
particles. Notice that if \i) is an eigenstate of S with a phase e 2lS , we have 

(i\T\B)* = e- 2lS ' (i\T\B) . (14) 

which implies equal rates for the charge conjugated decays and hence no CP asymmtery. Therefore, at least two 
different states with equal quantum numbers must exist which can be connected by strong rescattering. Also 

(i\T\B) = e u T t {i\T\B) = e l5 T* (15) 

The matrix elements Ti are assumed to be the "bare" decay amplitudes, calculated e.g. in factorization approxima- 
tion || and have no rescattering phases. This implies that these transition matrix elements between charge conjugated 
states are just the complex conjugated ones of each other. Eq. Jl5| ) is known as Watson's theorem 0. 
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The above considerations allow to formulate a condition for the non-vanishing CP asymmetry. It requires two 
different final states which undergo strong transitions into each other. The strong phase is then nothing but the 
occurence of the physical intermediate state \Bp) and arises when summing over the intermediate states. 

The final state rescatterings of high energy particles may be divided into 'soft' and 'hard' scattering. Soft scat- 
tering occurs primarily in the forward direction with limited transverse momentum, having a distribution which falls 
exponentially with a scale of order 0.5 GeV. Soft scattering might be best described by hadronic states. At higher 
transverse momentum one encounters the region of hard scattering, which falls only as a power of the transverse mo- 
mentum. Collisions involving hard scattering are interpreted as interactions between the quarks and gluons of QCD. 
Note that it is possible to generate FSI phases in nonleptonic B-decays into charmless final states in perturbative 
QCD [p|-p0|: there are two ways to reach a given final state, via the tree diagram b — ► uus, and via b — > ccs process, 
with subsequent final state rescattering of the two charmed quarks into two up quarks (penguin diagram). Since the 
energy release in b-decay is of the order mi, > 2m c , the rescattered c-quarks can go on-shell generating CP conserving 
phase and thus Aq p . 

For the soft FSI, the low energy effective theory of strong interactions can be used to reliably estimate FSI phase 
differences in the kaon system. In the D system final state rescattering has been studied assuming the dominance 
of intermediate resonances jn| . In the B system, where the density of the resonances available is large due to the 
increased energy, a different approach must be employed. One can use the fact that the 6— quark mass is large 
compared to the QCD scale. Then, the leading order behavior of soft final state phases in the — > oo limit can be 
investigated. This can be done by considering first the elastic channel, and demonstrating that elastic rescattering 
does not disappear in the limit of large tub ■ Since the unitarity of the 5-matrix requires that the inelastic channels are 
indeed the dominant contributors to soft rescattering, such contributions have to share a similar behavior in the heavy 
quark limit. The elastic channel is convenient because of the optical theorem which connects the forward (imaginary) 
invariant amplitude M. to the total cross section, 

1m M/->f(s, t = 0) = 2fc % /scr/^ a n ~ scr/^ a n , (16) 

where s is the squared center-of-mass energy and t is the squared momentum transfer. The asymptotic total cross 
sections are known experimentally to increase slowly with energy and can be parameterized by the form [[T2|fL3[ : 

\ 0.08 / \ -0.56 

S \ IS 



*<*>= X UJ +Y U • (17) 

where sq = 0(1) GeV is a typical hadronic scale. Thus, the imaginary part of the forward elastic scattering amplitude 
( |l6| ) increases asymptotically as s 108 . Considering only the imaginary part of the amplitude, and building in the known 
exponential fall-off of the elastic cross section in t (t < 0) H] by writing 



. 1.08 

S 



am M f ->f(s,t) ~ ifo I— j e bt , (18) 

one can calculate the contribution of the imaginary part of the elastic amplitude to the unitarity relation for a final 
state / = a + b with kinematics p' a + p' b = p a + pb and s = (p a + pb) 2 '- 



M B ^ f , (19) 



Vise Ms^f = 77 



167T sob \ Sq 

where t = (p a ~ p' a ) 2 — — s (l ~ cos0)/2, and s = 



1 .„ / 9 \ 0-08 
1 IPQ 
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On can refine the argument further, since the phenomenology of high energy scattering is well accounted for by the 
Regge theory |^,^5|. In the Regge model, scattering amplitudes are described by the exchanges of Regge trajectories 
(families of particles of differing spin) with the leading contribution given by the Pomeron exchange. Calculating the 
Pomeron contribution to the elastic final state rescattering in B — > 7777 one finds |16| 
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Vise MB-^nr|pomerou = -ieM B ^ww, £ ~ 0.21 . (20) 

It is important that the Pomeron-exchange amplitude is seen to be almost purely imaginary. However, of chief 
significance is the identified weak dependence of e on - the (m 2 ^) 008 factor in the numerator is attenuated by the 
\n.(m B / sq) dependence in the effective value of b. 

The analysis of the elastic channel suggests that, at high energies, FSI phases are mainly generated by inelastic 
effects. This conclusion immediately follows from the fact that the high energy cross section is mostly inelastic. This 
also follows from the fact that the Pomeron elastic amplitude is almost purely imaginary. Since the study of elastic 
rescattering has yielded a T-matrix element T a b^ a b = 2ie, i.e. S a b^ab = 1 — 2e, and since the constraint of unitarity 
of the 5-matrix implies that the off-diagonal elements are 0(y/e), with e approximately 0(m B ) in powers of tub and 
numerically e < 1, then the inelastic amplitude must also be 0(m B ) and of magnitude y/e > e. Similar conclusions 
follow from the consideration of the final state unitarity relations. 

The very presence of inelastic effects suggests a physical picture similar and complimentary to the "color trans- 
parency argument". This argument suggests that a "small" (compared to the typical hardonic size 1/Aqcd) color- 
singlet two-quark configuration does not interact with the soft gluonic environment. Thus, if the two-body decay 
is dominated by this particular quark configuration with all other quark configurations yielding multiparticle final 
states, one might expect that the effects of FSI are suppressed in the decays of ultra-heavy particles. However, this 
quark configuration is realized only on the edge of the available phase space. Therfore, in the limit mf, — * oo one ulti- 
mately encounters the situation where the quarks easily combine to form a multiparticle state which then undergoes 
rescattering into the two-body final state. Analysis of the final-state unitarity relations in their general form, 

ViscM B ^ h =\Y, M *^ T Lh > ( 21 ) 

k 

is complicated due to the many contributing intermediate states present at the B mass. However, it is possible to 
illustrate the systematics of inelastic scattering in a simple two-channel model. This example involves a two-body 
final state f 1 undergoing elastic scattering and a final state f 2 which represents 'everything else'. We assume that the 
elastic amplitude is purely imaginary. Thus, the scattering can be described in the one-parameter form 

f cos26< isin20\ (2is\n 2 9 sin20 \ , , 

b -\ism28 cos20 J ' ~ ^ sin20 2ism 2 ) ' [ ' 

where, from our elastic-rescattering calculation, we identify sin 2 9 = e. The unitarity relations become 

9 1 

Vise M B-*h = —i sin 6M b^j x + - sin 26M s_»/ 2 , 

Vise M B ->h = - sin 20M B ^ fl - i sin 2 9M B ^f 2 (23) 



Denoting M\ and Al!] to be the decay amplitudes in the limit 9^0, an exact solution to Eq. (^3() is given by 

M B ^h =cosOM 1 + ism6M% , M B->f 3 = cos 0M% + i sin 0M? • (24) 
In this example we see that the phase is given by the inelastic scattering with a result of order 



1m Mb^s r 



KeM B ^f ^ M\ ■ (25) 
Clearly, for physical B decay, we no longer have a simple one-parameter S matrix ,and, with many channels, cancella- 
tions or enhancements are possible for the sum of many contributions. However, the main feature of the above result 
is expected to remain - that inelastic channels cannot vanish because they are required to make the discontinuity real 
and that the phase is systematically of order y/e from these channels and thus does not vanish in the large tub limit. 
Moreover, it is possible to show that inelastic FSI can contribute to CP violating asymmetries at the leading order in 
tub [[l6| ■ Non-zero FSI phases have been recently observed by CLEO jl7| . 

We should note that radiative weak decays of B mesons are also affected by FSI. For instance, the extraction of Vtd 
from B — > p7 is hampered by uncertainties related to certain long distance effects from the on-shell hadron rescattering 
with subsequent conversion of one of the hadrons into the photon. This contribution could be sizeable |18| . 
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III. TESTING MODELS OF FINAL STATE INTERACTIONS 



(i) Model Independent Bounds on the FSI Corrections. In view of the large theoretical uncertainties involved in 
the calculation of the FSI contributions, it would be extremely useful to find a phenomenological method by which to 
bound the magnitude of the FSI contribution. The observation of a larger asymmetry would then be a signal for New 
Physics. Here the application of flavor SU (3) flavor symmetry, provides powerful methods to obtain a direct upper 
bound on the FSI contribution. 

The simplest example involves bounding FSI in B — > ttK decays using S ± — > K ± K transitions. The effective 
Hamiltonian for b — > d decays may be obtained from ([!]) by the substitution s — > d. In analogy with Eq. (|^) the 
amplitudes may be decomposed according to their dependence on CKM factors, giving 

A(B+ -> K+K°) = A cd - A ud e^e u , 

A(B- -» K-K°) = A cd - A ud e-^e lS . (26) 

Invariance under the SU(3) rotation exp(i^\j), i.e., interchange of s and d quark fields, implies equalities among 
operator matrix elements and amplitudes, 

(K-K°\Qf\B-) = (K\-\Qr\B-) ) g = u,c,; i = l,2 
(Jf"ii- |Qf|B-) = <Jf 7r-|Q?|S->, £ = 3,.., 10. (27) 
4je li = A+ e i5 +^(l + R ud ) , A cd = 4+^(1 + R cd ), 

*us Vcs 

where Rud and R c d parameterize SU (3) violation, which is typically of the order of 20 — 30%. Note that it is only 
an SU(2) subgroup of SU(3), namely [/-spin, which is required to derive these relations. Since the B~ carries U = 
and the transition operators Qj d and Qf carry U = |, it is only the U = h component of the K~K° final state which 
couples to the decay channel. As a result, an upper bound on e follows from the ratio 

R = BR(B+ -> K+ W) + BR(B- -> K_K°) 
K BR(B+ ^ K°tt+) + BR(B- ^ K°n-) 

After some algebra, we obtain for e and A^p 

e < Xy/R^(l + Re[R ud }) + \ 2 {R K + 1) cos 7 cos <5+ + C(A 3 , \ 2 Rud,cd) , 
\A% P \ < 2Xy/R K R(l + Re[R ud ]) + 2\ 2 VR{R k VT~ 7 R + R K + 1) + 0(\ 3 , \ 2 R ud , cd ) . (29) 

Using the experimental bound Rk < 1.9 one obtains e < 0.4 and Aq P < 0.6. More interesting constraints on e and 
^4qp could also be be obtained 

(izj Direct Observation. As we have shown above, the amplitude for the decay of a B- meson into some final state 
/ includes a direct contribution A(B — > f) and a sum over the contributions A(B — » % — > /), which corresponds to 
the weak decays of the S-meson into intermediate hadronic states i, followed by the strong scattering of i into /. 
The possibility of significant final state scattering effects becomes real when considering rare decays, for which the 
amplitude A(B — ► /) is suppressed compared to A(B — > i). In the ideal case, when the direct contribution is absent, 
one may be able to isolate the effect of FSI completely. While this situation might not be realized in the nature, rare 
weak decays offer tantalizing possibility of the direct observation of the effects of FSI. 

One of the possibilities involve dynamically suppressed decays which proceed via weak annihillation diagrams. It 
has been argued that final state interactions can modify the decay amplitudes, violating the expected hierarchy of 
amplitues. For example, it is expected that the amplitudes that do not involve spectator quarks (such as tree-level 
or penguin amplitudes) dominate over the diagrams involving spectator quarks (e.g. weak annihilation or weak 
rescattering amplitudes). In many cases, large amplitudes might contribute to the processes involving spectator 
quarks through the final state rescattering J20|]2l|] . It must be stressed that although the predictive power of available 
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methods is limited and most of the estimates are based on the two-body rescattering diagrams, some conclusions can 

still be reached. For instance, it is possible to show |2l| that the rescattering from the dominant channel leads to the 

suppression only of the order A ~ 0.2 compared to /b/tob ~ A 2 obtained from the naive quark diagram estimate. 

Another class of decays involve the so-called OZI-violating modes, i.e. those which cannot be realized in terms of 

quark diagrams without annihillation of at least one pair of the quarks. It includes modes like B d — > (fxj), D°(f> and 

— o 

J/ip<t>. For instance, the direct contribution to B d — > tficf) involves a suppressed space-like penguin diagram. However, 
the unitarity of the 5-matrix, S^S = 1, implies that this decay can also proceed via the OZI-allowed weak transition 
followed by final state rescattering into the final state under consideration. These two-step OZI violating processes 
were intensively studied in connection with certain low-energy processes ^2j,^3|. In i?-decays these OZI-allowed steps 
involve multiparticle intermediate states and might provide a source for significant violation of the OZI rule. For 
instance, the FSI contribution can proceed via B d -> tjMt/O -> 00, B d -> D* ?/') -> £>°0 and B d -► tp'r]^ -> J/tp<f>. 
The intermediate state might also include an additional set of pions. The weak decay into the intermediate state 
occurs at tree level, through the (uu + dd)/s/2 component of the rf'^ 1 wavefunction, whereas the strong scattering 
into the final state involves the ss component f24| . Hence the possibility of using these decay modes as direct probes 
of the FSI contributions to B decay amplitudes. It is however possible to show that there exist strong cancellations 
p3| , p5| among various two body intermediate channels. In the example of B d — > 00, the cancellation among r\ and r\' 
is almost complete, so the effect is of the second order in the SU (3)-breaking corrections 

Disc Mb^h> = 0{5 2 ,A 2 ,5A)f n F Q A, S = - /„ A = Fq — F , (30) 

with A ~ s Q0_1 e i7rQ °/ 2 /86. This implies that the OZI-suppressed decays provide an excellent probe of the multiparticle 
FSI. Given the very clear signature, these decay modes could be probed at the upcoming i?-factories. 
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